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ESTIMATION OF SEA SURFACE TEMPERATURE FROM 


REMOTE SENSING IN THE 11-13 WINDOW REGION 
C. Prabhakara, G. Dalu* and V. G. Kunde 
ABSTRACT 

The Nimbus 3 and 4 IRIS spectral data in the 11-13 p m water vapor win- 
dow region are analysed to determine the sea surface temperature (SST). The 
high spectral resolution data of IRIS are averaged over approximately 1 ptm wide 
intervals to simulate channels of a radiometer to measure the SST. In the 
present exploratory study three such channels in the 775-960 cm” 1 (12.9 - 
10.5 Mm) region are utilized to measure SST over cloud-free oceans. However 
two of these channels are sufficient in routine SST determination. The differ- 
ential absorption properties of water vapor in the two channels enable one to 
determine the water vapor absorption correction without detailed knowledge of 
the vertical profiles of temperature and water vapor. 

The feasibility of determining the SST is demonstrated globally with Nimbus 
3 data where cloud-free areas can be selected with the help of albedo data from 
the MRIR experiment on board the same satellite. The SST derived from this 
technique agrees with the measurements made by ships to about 1°C. 


*NRC-Research Associate. Permanent address: C.N.R.-Istituto di Fisica 
dell'Atmosfera - Rome, Italy. 
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ESTIMATION OF SEA SURFACE TEMPERATURE FROM 


REMOTE SENSING IN THE 11-13 ^m WINDOW REGION 
INTRODUCTION 

One of the techniques used to measure the sea surface temperature (SST) 
from satellites involves observations of the brightness temperature in the 8-13/./. m 
water vapor window region under cloud-free sky conditions. Allison and 
Kennedy (1967), Smith et al.(1970), Shenk and Solomonson (1972), Rao et al. 

(1972) and others have estimated the SST from remote sensing window measure- 
ments. The principal advantage with this infrared window is that it can be used 
both night and day. A second advantage is that low noise equivalent temperature 
may be attained in radiometric measurements as the terrestrial emission has its 
peak value in this spectral region. 

Several of the earlier attempts to estimate the SST from the window meas- 
urements have dealt elaborately with the problem of eliminating cloud contami- 
nated data. The present study does not offer any improvement in this respect. 

The problem of correcting the window data for water vapor absorption however, 
is better resolved in this study. 

The brightness temperature, measured from space over a cloud free ocean, 
in the window is in general lower than the SST. The water vapor absorption is 
mainly responsible for this effect. In the earlier studies the water vapor at- 
tenuation correction is derived from climatological data on water vapor and 
temperature in the atmosphere. Such climatological corrections can introduce 
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errors as large as ±2°C in the estimated SST. It is necessary therefore to ap- 
ply the water vapor correction corresponding to the local vertical distribution 
of temperature and water vapor to improve the accuracy of the estimated SST. 
Such soundings from a remote platform would require elaborate instrumentation 
and data reduction procedures. 

Anding and Kauth (1970), and Price (1973) have examined the usefulness of 
differential absorption in the window region to correct for the water vapor ab- 
sorption with soundings of temperature and water vapor. In this study a trans- 
mittance model for the water vapor window region is developed which can be 
used in a simple scheme that does not require temperature and humidity sound- 
ings of the atmosphere to estimate SST over cloud free oceans with the help of 
Nimbus 3 and 4 IRIS measurements. This scheme takes advantage of the differ- 
ential absorption properties of the water vapor in the spectral region 
775 - 960 cm" 1 . With this technique the SST is estimated with an accuracy of 
about 1°C from IRIS measurements. 

TRANSMISSION OF WATER VAPOR IN THE 11-13 M m REGION 

It is necessary to model the water vapor absorption in the three spectral 
regions 775 - 831, 831 -887, and 887 - 960 cm -1 which are used to simulate the 
broad band radiometer channels for the retrieval of SST. The 1000 - 1100 cm 1 
region has been avoided due to the strong ozone absorption. The effects of 
aerosols are not considered in this investigation. 
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The total water vapor transmission in the 11 jum window region may be 


expressed as 


T(V, T) = T e (v, T) T p (v, T) T^{V, T) (1) 

where r e and r p are associated with the water vapor continuum and Tj> with 
the water vapor lines. 

The laboratory measurements of Bignell (1970) and Burch (1970) show that 
the water vapor absorption in the 8 - 13 window region is primarily dependent 
on the water vapor partial pressure e. This absorption is referred to as 
"e-type” absorption, one of its characteristics being a strong negative tempera- 
ture dependence. Water dimer molecules may be responsible for this behavior. 
Window radiance measurements made by Platt (1972) with a radiometer on 
board an aircraft looking down at the sea surface and measurements made by 
Lee (1973) with a balloon borne radiometer, are consistent with the presence of 
the e-type absorption in the atmosphere. Nimbus 4 IRIS data reported by Kunde 
et al. (1974) also confirm this type of absorption. Grassl (1973) derived the 
magnitude of the aerosol and e-type absorption coefficients in the window region 
with the help of a ground based radiometer that measured the sky emission. 

The absorption coefficients deduced by Grassl are in close agreement with the 
laboratory data obtained by Bignell and Burch. Grassl also found that in tem- 
perate and cold climates, the aerosol absorption can be comparable to the e-type 
absorption. Wark (1972), Houghton and Lee (1972) and Prabhakara et al. (1972) 
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have emphasised the importance of "e" type absorption in determining the 
SST. 

The e-type absorption coefficients in the three channels considered here 
are obtained from the measurements of Burch (1970) made at a temperature of 
296°K. These absorption coefficients are assumed to decrease at a rate of 2% 
per 1°C (Bignell, 1970). That is 

k e(^. T) = k e (V, 296) [1 - 0.02(T - 296)]' ( 2 ) 

where k e is the e-type absorption coefficient (g -1 cm 2 ), 
v is the wave number (cm -1 ), 

T is the temperature (°K), 

The transmission is then calculated with the formula 

T e (v, T ) " Exp[-k e (V, T) ew] ( 3 ) 

where e is the water vapor pressure (atm), 

w is the path length of water vapor (g cm -2 ). 

The water vapor pressure in the atmosphere is assumed to decrease ex- 
ponentially with height with a constant scale height. This assumption implies 
that the mean vapor pressure e in a column of the atmosphere is proportional 
to the total precipitable amount of water w in that column. That is 
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e 


sw 


(4) 


where s is a constant. 

As a working hypothesis it is assumed that a mean vapor pressure of 3 mb 
in an atmospheric column of 1 cm 2 corresponds to 1 gram of precipitable water 
in that column, which means s = 3 (mb g _1 ). 

With these numerical values and with the help of equations (2) and (3), we 
have calculated the transmission r g for different path lengths at two tempera- 
tures 280° K and 300° K. 

The second component of the continuum absorption depends on the total 
pressure 


T p ( v < T) = Exp[-k p (y, T) wp/p Q ] 


(5) 


where 

kp is the absorption coefficient for the collision broadening in the wings of 
water vapor lines (g -1 cm 2 ), 

p is the average pressure of the layer of the atmosphere, weighted accord- 
ing to the distribution of water vapor, assumed to be 850 mb, 
p Q is the standard pressure of 1000 mb. 

Only upper limits have been inferred from the laboratory data for the value of 
k p . In this study we have adopted the values of k p used by Kunde and Maguire 
(1974). This component is only important for low water vapor concentrations 
(w < 1 g cm -2 ). 
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For the spectral lines. Goody (1954) has listed the line data of Benedict and 
Kaplan (1959) arranged in 20 cm' 1 intervals at three temperatures 220°, 260° 
and 300°K. These data are used to calculate the average transmission due 
to spectral lines, A statistical band model with non- over lapping lines, having 
an exponential line intensity distribution, has been assumed, 



where k^> = S/S is the mean absorption coefficient (g -1 cm 2 ), 

S is the mean line intensity (g 1 cm), 

S is the average line spacing (cm -1 ), 
a Q is the line half width at STP (cm - 1 ), 
p is 850 mb, 
p 0 is 1000 mb. 

The absorption coefficients used are listed in Table 1 for two temperatures, 
280° K and 300° K. The various transmission values derived in the three spec- 
tral intervals at the two different temperatures are shown in Table 2. The total 
transmission at the two temperatures is plotted in figs, la and lb as a function 
of path length for the three spectral channels. From these figures it is seen 
that the transmission function to a first approximation is linearly related to the 
path length up to about 4 grams of precipitable water at both temperatures. 

This approximate linearity of the transmission is also obtained when we assume 
that the line intensity distribution has a form proposed by Godson (1954). As 


6 



Table 1 


Water vapor absorption coefficients 


Channel 

No. 

Av (cm” 1 ) 

k p (g- 

1 cm 2 ) 

k e (g~ 

1 cm 2 ) 

k i <s" 

cm 2 ) 

a 0 /S 

280° K 

300° K 

280° K 

300° K 

280° K 

300° K 

STP 

1 

775-831 

0.035 

0.040 

19.46 

13.56 

0.333 

0.497 

0.015 

2 

831-887 

0.017 

0.020 

14.52 

10.12 

0.098 

0.157 

0.018 

3 

887-960 

0.009 

0.010 

11.59 

8.08 

0.047 

0.074 

0.014 
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Table 2 A 


Water Vapor Transmissivities for 280° K 


Channel 




Precipitable water (g cm" 

2 ) 



0.5 

1 

2 

3 

4 

6 

8 

( 

1 

775- 

831 

0.985 

0.970 

0.942 

0,914 

0.887 

0.835 

0.787 

Tp 1 

2 

831- 

887 

0.993 

0.986 

0.971 

0.957 

0.944 

0.917 

0.890 

1 

3 

887-960 

0.996 

0.992 

0.985 

0.977 

0.970 

0.955 

0.940 


r 

1 

775- 

-831 

0.985 

0.943 

0.792 

0.591 

0.393 

0.122 

0.024 

T e < 

2 

831- 

-887 

0.989 

0.957 

0.840 

0.676 

0.498 

0.208 

0.062 



887-960 

0.991 

0.966 

0.870 

0.731 

0.573 

0.286 

0.108 


r 

1 

775- 

-831 

0.920 

0.880 

0.825 

0.783 

0.748 

0.690 

0.641 

T l < 

2 

831 

-887 

0.963 

0.939 

0.905 

0.878 

0.856 

0.820 

0.790 


3 

v_ 

887 

-960 

0.981 

0.967 

0.947 

0.931 

0.917 

0.895 

0.877 


1 

775 

-831 

0.893 

0.805 

0.615 

0.423 

0.261 

0.070 

0.012 

r < 


831 

-887 

0.946 

0.886 

0.738 

0.568 

0.402 

0.156 

0.044 



CO 

00 

-960 

0.968 

0.927 

0.812 

0.665 

0.510 

0.244 

0.089 
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Table 2B 


Water Vapor Transmissivities for 300° K 


Channel 

tu 


Precipitable water (g cm" 

2 ) 


0.5 

1 

2 

3 

4 

6 

8 


r 

1 

775-831 

0.983 

0.967 

0.934 

0.903 

0.873 

0.815 

0.762 

T < 
P 1 

2 

831-887 

0.992 

0.983 

0.967 

0.950 

0.934 

0.903 

0.873 


3 

v_ 

887-960 

0.996 

0.992 

0.983 

0.975 

0.967 

0.950 

0.934 



775-831 

0.990 

0.960 

0.850 

0.693 

0.521 

0.231 

0.074 

T e 

2 

831-887 

0.992 

0.970 

0.886 

0.761 

0.615 

0.335 

0.143 


3 

887-960 

0.994 

0.976 

0,908 

0.804 

0.678 

0.418 

0.212 


fl 

775-831 

0.899 

0.851 

0.784 i 

0.733 

0.691 

0.620 

0.560 

T i < 

2 

831-887 

0.948 

0.917 

0.874 

0.841 

0.814 

0.769 

0.731 


3 

V. 

887-960 

0.972 

0.954 

0.928 j 

0.909 

0.892 

0.865 

0.842 


1 

775-831 

0.875 

0.790 

0.622 

0,459 

0.314 

0.117 

0.032 

T < 

2 

831-887 

0.933 

0.874 

0.749 

0.608 

0.468 

0.233 

0.091 


3 

887-960 

0.962 

0.924 

0.828 

0.713 

0.585 

0.343 

0.167 
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w (g cm-2) 


Figure 1. Water vapor transmission in the three channels; 
1. 775-831 cm' 1 , 2. 831-887 cm' 1 and 3. 887-960 cm r ’ ; as £ 
function of precipitable water content in the atmosphere for 
T = 280° K (a) and T = 300° K (b). 





these temperatures and path lengths cover a wide variety of atmospheric condi- 
tions from tropics to high latitudes, it may be stated that this approximate linear 
law is applicable over most of the globe. Using the linear property of the trans- 
mission function shown in figs, la and lb, the following equation is developed 

T(v, w, T) “ 1 - K(V) w(l - e(T - 296^ ( 7 ) 

where K(v ) is an effective absorption coefficient corresponding to the total 
water vapor absorption and the coefficient e accounts for the temperature 
dependence of K(i'). 

The transmission functions derived here are dependent on the assumed 
water vapor distribution in the atmosphere, on the statistical band model, and 
on the absorption coefficients and k p . For this reason the magnitude of K(> ) 
can change with the choice of the model and the parameters used in the calcu- 
lations. However it is found that for a realistic combination of parameters, the 
transmission function r (v) varies in a nearly linear fashion with respect to 
path length w. This property is exploited in the present study. 

THEORY 

In a cloud-free, non scattering atmosphere under local thermodynamic 
equilibrium, the radiative transfer equation for the upwelling intensity I (v ) may 
be written as 
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l(v) = B(v, T s ) T(v, p 0 ) + I B[v, T(P)] dr(v, p) (8) 

P Q > 

where p 0 and p are surface pressure and pressure at any height (mb) respectively, 
v is the wave number (cm -1 ), 

T is the temperature (K), 

T s is the surface temperature (K), 

B is the Planck intensity (erg cm -1 ster -1 s _1 ), 

r is the transmission from any pressure level p to the top of the atmosphere. 
The surface emissivity in equation (8) is assumed to be unity; 

The above equation may be simplified as 

I(v) = B(v, T s ) r(y t p Q ) + B (v) [1 - r(v, p Q )] (9) 

where B(v) is tlie weighted mean Planck emission of the atmosphere 

l 

dr (10) 

V 

Substituting equation (7) into equation (9), we get 

I O') = B(v, T s ) - AK(v) [B(v, T s ) - B (u)} < n > 

where A = w [1 - e(T(p) - 296)] is a constant for a given atmosphere, independent 
of wave number. 


B(v) = 


■j: 

J 7 (V , 


B [v, T(p)] dr 




/I 
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The Planck function B can be expanded about the surface temperature T s 
with only the linear terms retained to yield 


3B(v, T ) 

B(v, T) * B(v, T s ) +■ (T - T s ) 


(12) 


This approximation holds good over a small range of temperatures and a narrow 
wave number span. With this approximation equation (11) becomes 


T(V) = T, - AKO) [T, - TO)] (13) 

where TO) is the brightness temperature and TO) is the equivalent radiative 
temperature of the atmosphere. 

Eq ua tion (13) indicates a linear relationship between the observed brightness 
temperature and the absorption coefficient KO), provided that T(y) is not 
strongly dependent on v. With the help of the radiative transfer equation (8), 
model calculations based on typical atmospheric temperature and humidity pro- 
files were made to examine the amount by which T(v) changes with v. It is 
found from these model calculations, that within the 775-960 cm -1 region, TO) 
changes with v by less than 1°C. The magnitude of (T - T (v )) is found to be 
about 10°C. In view of this, the change in TO) is neglected in equation (13) to 
arrive at an approximate linear relationship between TO) and KO) as follows 

TO) = T, - ySKO) ( 14 > 

where P = A [T - TO)] is a constant for a given atmosphere. 
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Now if we know the brightness temperature T(v) in at least two spectral 
regions, which have appreciable difference in K(v), equation (14) can be used 
to estimate the sea surface temperature T s . This formulation does not require 
explicitly either the temperature and humidity distribution in the atmosphere, 
or the temperature and pressure dependence of the absorption coefficients. 

RESULTS AND DISCUSSION 

The Nimbus 4 IRIS had a-^ 95 km field of view, 2.8 cm' 1 spectral resolution 
and noise equivalent radiance of about 0.5 erg cm -1 ster' 1 s -1 (Hanel et al„ 

1972). These high spectral resolution data reveal the influence of numerous 
water vapor absorption lines, which are strong in the tropics and decrease in 
strength at high latitudes. In addition these spectra indicate that the brightness 
temperature increases steadily from about 775 cm' 1 to 960 cm" 1 . This is il- 
lustrated in Fig. 2 where a few spectra selected over cloud-free ocean are 
shown. The cloud-free condition is inferred from the Nimbus 4 Image Dissector 
Camera System (IDCS) pictures (Sabatini, 1970). This procedure is somewhat 
subjective and is not practical for collecting a large sample of cloud-free spectra. A 
small sample of 8 clear sky spectra was selected where sea surface tempera- 
ture was available from ship measurements. The mean brightness temperature 
corresponding to each one of the spectral channels 775-831, 831-887 and 
887-960 cm" 1 is listed in Table 3. From this sample of IRIS data and the 
corresponding observed sea surface temperature, one can determine the relative 
values of the three absorption coefficients K(v) from equation (14) for each of 
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Table 3 


Nimbus 4 IRIS window brightness temperature, 
and the SST from IRIS and ships. 


No. 

Lat. 

Long. 

T 

l 

T 

2 

T 

3 

T 

IRIS 

T 

SHIP 

1 

41.0°N 

62.5°W 

272.9 

275.2 

276.8 

281.2 

281.0 

2 

36 .6°N 

60.9°W 

284.5 

286.9 

287.9 

292.0 

1 

3 

11.0°N 

53.5°W 

287.8 

291.7 

293.4 

300.1 


4 

35.6 D N 

33.8°W 

281.2 

283.8 

285.0 

289.6 

291.0 

5 

44.6°N 

37.1 °W 

280.9 

283.1 

284.0 

287.7 

287.0 

6 

19.7°N 

82.6°W 

291.2 

294.2 

295.6 

300.7 

301.5 

7 

15.1°N 

144. 7°W 

292.7 

294.9 

296.1 

300.1 


8 

39.4°N 

69.1°W 

286.9 

290.6 

291.8 

298.0 

296.2 


Table 4 

Relative absorption coefficient (g _1 cm 2 ) for the 
three window channels 


Channel 

1 

2 

3 


775-831 

831-887 

887-960 

K(v) 

0.191 

0.131 
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the three spectral intervals, A first guess for these coefficients is derived 
from the synthetic transmission function illustrated in Figure la and these 
values are then adjusted to give a best fit to the data of Table 3. These values 
of K(v) are listed in Table 4, In Fig, 3 a plot of the measured brightness tem- 
peratures versus K(i>) for the 8 cases is presented to show the degree of 
linearity in the observed data. 

The advantage of using the method presented here to determine the SST is 
that it is inherently capable of taking into account the effects of water vapor and 
temperature profiles satisfactorily. This capability is particularly necessary 
in measuring the SST with meaningful accuracy over the tropics where the water 
vapor absorption correction is large and quite variable in time and space. To 
demonstrate this point a set of 6 cloud free Nimbus IRIS spectra from one orbit, 
from 8°N to 14°N along 55°E longitude on May 10, 1970, over Arabian Sea, are 
analysed. The window brightness temperature in the three channels for each of 
these 6 spectra is plotted in Fig. 4, as a function of the absorption coefficient 
K( v). The intercept on the ordinate gives the value of the SST, It may be noticed 
that the IRIS window brightness temperature in the most transparent channel 
(887-960 cm -1 ) increases by about 3°K going northward from 8°N to 14°N. One 
would erroneously infer from water vapor correction based on climatological 
atmospheric temperature and humidity data an increase, in the SST, of about 
3°K from 8°N to 14°N. However, with the present correction scheme the esti- 
mated SST does not increase towards the north; instead it shows a small 
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Figure 3. IRIS window brightness temperature versus relative absorption coefficient. 
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Figure 4. Window brightness temperature variation over the Arabian Sea be 
tween 8° N to 14° N along 55° E longitude on May 10, 1970, from Nimbus 4. 




decrease of about 1° in agreement with the seasonal pattern (U.S. Naval Ocean- 
ographic Office, 1967) for that oceanic region. 

The Nimbus 3 IRIS had a~150 km field of view, 5 cm” 1 spectral resolution 
and noise equivalent radiance of about 1 erg cm” 1 ster -1 s -1 (Hanel and Con- 
rath, 1969). These data, although slightly inferior to the Nimbus 4 IRIS meas- 
urements, can serve as good check on the validity of the method and the magni- 
tude of the coefficients K(v) because we can objectively select cloud free data 
in the following fashion. Measurements from the Nimbus 3 satellite, of a scanning 
Medium Resolution Infrared Radiometer (MRIR), with a field of view of about 50 km 
at nadir, included the reflected solar radiation in the visible and the infrared radi- 
ation in the 11 Mm region (Sabatini, 1969). Clouds can reflect solar radiation in 
significant amounts while the ocean surface cannot. Thus with the help of the 
MRIR reflected solar radiation measurements at the nadir, it is possible to 
objectively select cloud-free IRIS data. During a period of 10 days in April 1969, 
a sample of 106 IRIS measurements over the oceans was extracted when corre- 
sponding ship measurements of the sea surface temperatures were avail- 
able within ±1° of latitude and longitude. These data covered the globe from 
60° S to 60° N. The SST derived from these IRIS data are compared with the 
ship data in Fig. 5. We see that over a wide dynamic temperature range of 
about 40° F to 85° F, the two sets of data agree very well. The R.M.S. difference 
between the two sets of data is about 2.5° F or 1.3° C. However the SST meas- 
ured by ships has an error of about 1° C (Saur, 1963) and so the IRIS-SST 
measurements should have an accuracy better than 1.3°C. 
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In order to examine if there is any systematic bias in the IRIS-SST with 
respect to ship measurements, a linear regression analysis is performed to 
obtain a line of best fit as shown in Fig. 5. One can notice that this regression 
line does not coincide with the line of perfect agreement which would be inclined 
at 45° to the axes. It appears that the use of K(^) given in Table 4 leads to a 
slight underestimation of the SST over cold water of high latitudes and an over- 
estimation of the SST in the tropics. A further refinement in the magnitude of 
K(i') eliminates this systematic bias and improves the accuracy of determining 
SST to about 1° C. 

CONCLUSION 

The basic result of this study is that over cloud free oceans we can make 
use of the differential water vapor absorption properties of the 11 p m window 
region to correct for the atmospheric attenuation. This correction scheme in- 
herently takes into account the temperature and humidity structure prevailing at 
the location. Although three measurements in the window region 775-960 cm -1 
were used to demonstrate this method, in principle two measurements are 
sufficient. The two window measurements may preferably be in the approximate 
regions 775-831 and 887-960 cm" 1 to obtain the maximum effect of differential 
absorption. 

The spatial and temporal changes of sea surface temperature in the tropics 
are generally small (about 2 ° C). The water vapor absorption correction over 
the tropics could range from about 4° to 8° C. Hence the technique of estimating 
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the SST proposed in this study would be particularly valuable over the tropical 
oceans . 

As the radiometric errors are amplified in the derived results, it would be 
desirable to have a radiometer capable of measuring brightness temperature 
accurate to 0.1° K in order to estimate the SST with an accuracy of about 
0.5° K. 

The scheme presented here can be adopted to any existing sea surface 
temperature measurement program, based on one window channel, by adding 
one more channel in the window. 
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